AN IZERGIN KOREPIN-TYPE IDENTITY 
FOR THE 8VSOS MODEL, WITH APPLICATIONS TO 
QQ ■ ALTERNATING SIGN MATRICES 

o 
o 

^ ■ HJALMAR ROSENGREN 

>>: 

Abstract. We obtain a new expression for the partition function of the 8VSOS model with 
domain wall boundary conditions, which we consider to be the natural extension of the Izergin— 
Korepin formula for the six-vertex model. As applications, we find dynamical (in the sense of 
' the dynamical Yang-Baxter equation) generalizations of the enumeration and 2-enumeration 

of alternating sign matrices. The dynamical enumeration has a nice interpretation in terms 
of three-colourings of the square lattice. 

o 
u 

i 1. Introduction 

An alternating sign matrix is a square matrix with entries 0, —1 and 1, such 
that the non-zero entries in each row and column form an alternating sequence of 
^ I the form 

o>: 1,-1,1,-1 -1,1. 

CN ■ Mills, Robbins and Rumsey |MRR] conjectured that the number of n x n alter- 

nating sign matrices equals 

O' , l!4!7!---(3n-2)! 

00 : An ^ ' 



O! n\{n + 1)\ ■ ■ ■ {2n - 

This was proved thirteen years later by Zeilberger fZej . Kuperberg ^Klj found 
^ . a simpler proof based on the six-vertex model on a square with domain wall 

boundary conditions. This is a lattice model of statistical mechanics, whose states 
can be identified with alternating sign matrices. Each state carries a weight, in 
general depending on 2?7, -f- 1 parameters q,Xi, . . . ,Xn,yi, ■ ■ ■ ,yn- The partition 
function for the model is the sum of the weight of all states. By the Izergin- 
Korepin identity [H lICKj , it can be expressed in terms of the determinant 



det 



1 



i<i,j<n \ {xi - qyj){xi - q ^yj] 

Kuperberg observed that when q = e^^^^^ and Xi = yi = 1 for all i, the weight 
of each state can be normalized to 1, so the partition function is equal to An- The 
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Izergin-Korepin identity then gives 



A„ = 3("2') lim TT 

xi,...,x„^l -'--L 



yi,...,yn^l l<«<i<" 



{xi - Xj){yi - Vj) i<i,j<n \xf - 




) 



Although computing the hmit is not trivial, Kuperberg could do it by elementary 
means. 

The eight-vertex model is a generalization of the six-vertex model, where the 
weights may be taken as elliptic functions of the parameters. In his solution of the 
eight-vertex model, Baxter |B2] introduced a different generalization of the six- 
vertex model, the 8VS0S (eight-vertex-solid-on-solid) model. Actually, the "8" 
is somewhat misleading, since the model obeys the ice rule and thus admits only 
six local configurations. In particular, imposing domain wall boundary conditions, 
states can be identified with alternating sign matrices. Compared to the six-vertex 
model, its main distinguishing feature is the presence of a "dynamical" parameter, 
associated to the faces of the underlying lattice. 

The purpose of the present paper is to study generalizations of the Izergin- 
Korepin identity and of Kuperberg's specialization, when the six-vertex model is 
replaced by the 8VS0S model. One motivation is to understand the significance 
of the dynamical parameter from a combinatorial viewpoint. We also hope that 
our results may be useful for studying the thermodynamic limit of the partition 
function, similarly as the Izergin-Korepin identity is used in \KZ\ [Zi] . 

The plan of the paper is as follows. §2] contains preliminaries on theta functions, 
and in ^we recall the 8VS0S model and obtain some elementary properties of its 
partition function. In §l]we show that the partition function can be identified with 
a special case of the elliptic weight functions of Tarasov and Varchenko |TVj . We 
then give our main result. Theorem 15.11 which expresses the partition function 
as a sum of 2^ determinants. We argue that this is a natural extension of the 
Izergin-Korepin identity. The rest of the paper is concerned with the case when q 
is a root of unity. In ^ we show that if = 1 , our generalized Izergin-Korepin 
identity can be reduced to a sum of — 1 determinants. After the preliminary 
§71 we consider the analogue of Kuperberg's specialization for the trigonometric 
8VS0S model in ^ Curiously, the corresponding extension of the alternating 
sign matrix theorem. Corollary 18.41 contains not only the numbers An, but also 
the numbers 



which enumerate cyclically symmetric plane partitions in a cube of size n (Xj. 
This dynamical enumeration is best understood in terms of three-colourings of the 
square lattice. In particular, it allows us to compute exactly the probability that 
a random square from a random three-colouring, satisfying domain wall boundary 
conditions, has any given colour, see Corollary 18.51 Finally, in ^ we study the 
more elementary case of 2-enumeration. 



n 



(3j-l)(3j-3)! 
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During the completion of our work, several closely related papers have appeared. 
Pakuliak, Rubtsov and Silantyev |PRSj studied the 8VS0S partition function us- 
ing algebraic techniques, and independently obtained Theorem 14. 1[ Foda, Wheeler 
and Zuparic [FWZj studied some other elliptic models, where the partition func- 



tion can be explicitly factored. Proposition 19.11 can also be obtained from their 
results. Finally, Razumov and Stroganov |RSj studied the partition function of 
the three-colouring model using an elliptic function parametrization. Trying to 
understand the relation to their work greatly improved the exposition in §SJ 

Acknowledgement: I would like to thank Professors Michio Jimbo, Alexander 
Razumov, Yuri Stroganov, Vitaly Tarasov and Ole Warnaar for vital comments 
at various stages of the work. 

2. Theta functions 

Throughout, r and t] will be fixed parameters such that Im(r) > 0, t] ^ Z + rZ. 
We will write p = e^'^^'^ and q = e^^^^. By we always mean e^^^^^. 
We will use the notation 

oo 

[x] = - - p^+ig-"). 

j=0 

Up to a multiplicative constant, [x] equals the Jacobi theta function 6i{rix\T) 
[WW]. We sometimes write for short 

[Xi, . . . ,X„] = [Xi] ■ ■ ■ [Xn]. 

The function x ^ [x] is odd, entire, and satisfies 

[x+u,x — u,y + v,y — v] — [x + v,x — v,y+u,y — u] = [x+y,x — y,u+v,u-~v]. (2.1) 

In fact, up to an elementary multiplier, the only such function is the Jacobi theta 
function, together with the degenerate cases [x] = sin(7r?7x) and [x] = x, cf. |WWt 
p. 461]. We find it helpful to think of [ two-parameter deformation of the 

number x. 

In the case = 1, we find it more convenient to use the notation 

oo 

e{x) = e{x;p) = Y[{l-p^x){l-p'+^/x), 

j=0 

9{xi, ...,Xn) = 9{xi, . . .,Xn;p) = 9{xi;p) ■ ■■9{xn;p), 

so that 

[a;]=g--/2^(g-p). 
The following terminology will be useful. 

Definition 2.1. Fixing r and rj, we say that f is a theta function of order n and 
norm t if there exist constants ai, . . . , a„ and C with ai + ■ ■ ■ + an = t, such that 

f{x) = C[x-a^]---[x-a^]. (2.2) 
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Equivalently, f is an entire function such that 

fix + 1/r]) = fix + t/t]) = (_i)«e2-^''(*-"^)— ^™/(a;). (2.3) 

The equivalence of these two properties is classical |Wet p. 45]. More generally, 
any function of the form 



m+nj 



^^^^ ^ [x-b['^]---[x-b^'^ 



m 



where 

for each j, satisfies the quasi-periodicity (12. 3p . If / is entire (that is, the singu- 
larities at X = b['^'^ are all removable), it can then be factored as in (12. 2p . Unless 
/ is identically zero, the zero set is then + Z?7~^ + Ztt]~^, 1 < i < n, where 
ai + ■ ■ ■ + ttn = t. Thus, to prove that / vanishes identically, it suffices to find 
n independent zeroes. This gives a powerful (and classical) method for proving 
theta function identities, which we are going to apply repeatedly. 
Finally, we recall Frobenius' determinant evaluation jF?] 



Xi 



y, + i] 



det 

l<hj<n \ [Xi - yj\ 



-lp[tr-^[\x\-\y\+t] 111 



(2.4) 



or equivalently 

'eitXi/yj] 



det 



i<i,j<n \ 6ixi/yj] 



i-lp9itr-'Y9itx/Y)u,<^<,<n^,yA^^/^vy^/yJJ 



(2.5) 



Illj=lyj0i^^/y3) 

Here and throughout, we write 

\x\ = Xi + ■ ■ ■ + Xn, X = Xi---Xn- (2.6) 

3. The 8VS0S model 

We will study the 8VS0S model on a square with domain wall boundary condi- 
tions. There are several ways to describe this type of models, see e.g. [P]. We find 
it convenient to use height matrices. Fixing a non-negative integer n, by a state 
we mean an (ra + 1) x (ra + 1) matrix, such that any two horizontally or vertically 
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adjacent entries differ by 1, and sucfi tliat tlie boundary entries are specified as 

1 2 



1 
2 



n \ 

n — 1 
n-2 



(3.1) 



\n n — 1 n — 2 ■■■ / 

As an example, wlien n = 2 tliere are two states: 





(3.2) 



By a block, we mean a (2 x 2)-block of adjacent matrix entries in a state. Tlie 
blocks can be viewed as entries of an n x n matrix. Replacing each block ( ^ ^ ) 
by {b + c — a — d)/2 gives a bijection between states and alternating sign matrices 
of size n X n. For instance, the states (13. 2p correspond to the alternating sign 
matrices 

We will consider Boltzmann weights labelled by a,b,c,d G {±1} = sat- 
isfying the "ice rule" a + b = c + d. For any such labels, let there be given a 
meromorphic function i?"^ of two complex variables. We also fix 2r;, + 1 generic 



complex parameters A,a;i, 



If a block 



c d . 



has coordinates (i, j), chosen with 1 < z,j < n in the standard 



way, that block is said to have local weight 

j-,b—a,d~b / \ I \ 

This describes a generalized ice model, where "generalized" refers to the "dynam- 
ical" or "face" parameter A, which is absent in the six-vertex model. The weight 
of a state is defined as the product of all local weights, and the partition function 
as 



Znix; y; A) 



E 

states 



weight (state). 



(3.3) 



As an example, from (13. 2p we see that 

Z2{x;y; A) = Rt+{X,Xi - yi)i?+l(A + - ya) 

X RlX{\ + l,X2-yi)Rt+{\x2-y2) 
+ RXX{\ xi - yi)Rt+{X + 1, si - y2) 

X (A + 1, X2 - yi)RZZ (A + 2, - 1/2) • 

One is particularly interested in Boltzmann weights satisfying the quantum dy- 
namical Yang-Baxter equation (or star-triangle relation), which can be described 
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as follows [F]. Let V = Ve+ © Ve^ be a two-dimensional complex vector space, 
and introduce the operators R{\,u) G End(V" (S> V) by 

c+d=a+b 

Then, 

i?^2(A + /^^ Ml - M2)i?^^(A, Ui - U3)i?^^(A + M2 - Ms) 

= i?'=^(A, U2 - U3)i?'^(A + h\ ui - U3)R'^{X, ui - U2). (3.4) 

This should be understood as an identity for meromorphic functions of A, ui, U2, U3 
with values in End(V"®^), with notation as explained by the example 

i?^^(A + h^, ui - U2){ea ® eb® Cc) = R{\ + c, ui - U2){ea ® Cb) ® Cc- 

Part of the interest in this case comes from the following fundamental fact, which 
follows from the discussion in |B3l §9.6]. 

Proposition 3.1 (Baxter). If the Boltzmann weights satisfy (13.41) . then the par- 
tition function Zn{x; y; A) is a symmetric function of x and y. 

Adopting the normalization of [D], the 8VS0S model is given by the following 
solution of ([33D: 

[u + 1] 



Rlti^,u) = RZZ{X,u) 



[1] 



R+-i\u) - ^Jpp, ^-+(A,«) - ^ipp: 

i?;l(A,n) = l^, te(A,«)-[^-"] 



[A] ' ' ^ [A] ■ 

From now on, we restrict our attention to this model. 

Lemma 3.2. The partition function is a theta function of each Xi of order n and 
norm \y\ + A, and of each yi of order n and norm \x\ — A. 



Proof. By Proposition 13. ![ it suffices to consider the case i = 1. Thus, we consider 
/ as a function of xi] the case of yi is treated similarly. It is well-known and easy 
to see that for each state there exists a fc, with 1 < /c < n, such that the second 
row of the height matrix is 

123 ■■■ kk-lkk + l ■■■ n - I. 

The Xi-dependent part of the partition function is then 

[xi-yi + l]--- [xi - + 1][A + A; - 1 - a;i + yk][xi - yk+i] ■■■[xi- yj, (3.5) 

which is a theta function of the desired form. □ 



IZERGIN-KOREPIN IDENTITY FOR 8VSOS MODEL 



7 



Lemma 3.3. The partition function satisfies 

_ [A + n] Yll=2[y^ -Vk- l][xk - yi] 



Zn{x; y; A) 



Znix; y; A) 



xi+i=yi [A + r2-l][l]2("-i) 

X Z„_i(x2,...,a;„;?/2,---,2/n;A), 

^ nL2bi -yk + -1/1 + 1] 
X 2'„_i(x2, . . . ,a;„; ?/2, • • • ,2/n; A + 1) 



By Proposition 13.11 there are similar identities for any specialization Xi = yj 

and Xi = yj + 1, 1 < i,j < n. 



Proof. We first consider the case xi + 1 = yi. We observe that if 7^ 1 , then fl3.5p 
vanishes. Thus, only states with k = 1 contribute to the partition function. This 
fixes also the second column, all such states having the form 



1 
2 
3 



1 

1 
2 



2 
1 



3 
2 



n \ 
n — 1 
n-2 
n — 3 

J 



\n n — 1 n — 2 n — 3 ■ ■ ■ 
It follows that the partition function factors as 

n 

Rt+{X, -1) n ^^=(^ + ^ - 1' ^1 - yk)RzX{X + k-l,xk-yi) 

X Zn-l{x2, . . . ,Xn;y2, ■ ■ ■ , yn] A), 



k=2 



which simplifies to the desired expression. 

To prove the second identity it is better to let xi = yn, which is equivalent in 
view of Lemma 13. 3[ One may then apply a similar argument as before. □ 



4. Relation to elliptic weight functions 

In [TV], Tarasov and Varchenko introduced elliptic weight functions, which 
have played a fundamental role for constructing solutions to the gKZ and gKZB 
equations, see further |FTVH[FTV2] . The following result shows that the partition 
function for the 8VS0S model is an elliptic weight function. Indeed, if we let £ = n 
and ^1 = ■ ■ ■ = = y/v in |TVt Eq. (2.20)], it is straight-forward to identify the 
two expressions. 
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Theorem 4.1. The partition function can be represented as 



Zn{x; y; A) 



Ulj=i[yj - 




(4.1) 



Proof. Consider the two sides of (14. ip as functions of xi. By Lemma [3. 2[ the left- 
hand side is a theta function of order n and norm \y \ + A, and it is straight-forward 
to check that the same is true for the right. Thus, as discussed in ^ it suffices 
to verify fl4.ll) for Xi = yj — 1, 1 < j < n. Since both sides are symmetric in y, it 
is in fact enough to take xi = yi — 1. On the left-hand side, we may then apply 
Lemma 13.31 On the right, only terms with a{l) = 1 are non-zero, so it can be 
viewed as a sum over Sn-i- In this way, fl4.ip is reduced to the same identity with 
n replaced by n — 1, and is thus proved by induction on n. □ 

The same result was independently obtained by Pakuliak et al. |PRSj . 



The following identity is our main result. Originally, we derived it from The- 
orem 14.11 by a complicated argument. However, as soon as one has guessed the 
formula, it is easy to prove directly. 

Theorem 5.1. For generic 7, the partition function can be represented as 



5. An extension of the Izergin-Korepin identity 



Zn{x;y; A) 




Ui<i<j<n[^i - Xj][yi - Vj] 




[xf - Vj + 7] 




where 




(5.2) 



Before we prove Theorem 15. H we make some remarks. 
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Remark 5.2. The determinants in (15.11) are evaluated by (12.41) . This leads to the 
equivalent identity 

7 ( \\ [A + n] 

Z^[x] y; A) - 



[l]"'[7][|x| - |?/| + A + 7 + n] 



^ , MSI [A + 7 + ^ - l>^|] [\A - + 7 + \S\] 

Y [^ + n-\S\] 

X n H^^u(Ui-^-yM-^ + '-y^])' (5-3) 

which expresses Z„ as a sum of 2" explicitly factored terms. This is much better 
than the terms in (13. 3p or the n\ terms in (14. ip . In the case of the six- vertex 
model, the corresponding identity (without the freedom of choosing 7) is discussed 
by Warnaar [Wl Eq. (3.3)]. 

Remark 5.3. Theorem 15.11 can be viewed as a generalization of the Izergin- 
Korepin identity. To see this, consider the degenerate case when [x] = sin(7rr7x), 
with Im(?7) > 0, and A — 00. Up to normalization, this limit corresponds to the 
six- vertex model. Since [A + a]/[A + 6] ^'^'^vib-a) ^ obtain in the limit 

^ , , (_l)(;j)eWN-l.l) ^^=Jx. - y,][x, + 1 - y,] 

Zn{x;y; 00) = - 



[1]" [7]" Ul<^<,<n[^^-Xj][yi-yj] 

V (-1)1^1 det (tl^yilJi] 

ff^.X l<M<nV [Xf-yj] J 



l<i<j<nl 
X 



5C{l,...,n} 

By linearity of the determinant and a trigonometric identity, the sum in S can be 
written as the single determinant 



det 



Xi 



Vj + 7] [xi + 1 - yj + 7] 



i<i,j<n \ [xi - yj] [xj + 1 - yj] 



i<i,j<n \ [xi - yj\ [x, + l- yjl 



Thus, 7 cancels and we obtain the Izergin-Korepin identity in the form 
(_l)(^)eWkl-l.l) n-,=ik. - yj][xi + 1 - yj] 



Znix;y;oo) 



Ul<^<,<n{^^-XJ][y^-yJ] 



^ i<l?<n ( [xi - yj] [xi + 1 - y-j] 
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Proof of Theorem \5.1[ When n = 1, the result can be simphfied as 

[\ + j + l][x-y + j][\][x-y + l] - [\ + j][x-y + -f + l][\ + l][x-y] 

= [l][-f][y - X + \][x - y + \ + + I], 

which is an instance of (12. ip . 

We now proceed by induction on n. Let /l and ffj denote the left-hand and 
right-hand sides of (15.11) multiplied by [\x\ — \y\ + X+'y+n], and viewed as functions 
of yi. It is easy to verify that both fi and f^ are theta functions of order n + 1 
and norm 2\x\ + j + n — y2 — ■ ■ ■ — yn- That the singularities at [yi — yj] = are 
removable follows from a symmetry argument, and is also apparent from (15.31) . 
Thus, it suffices to verify (15.11) for n + 1 independent values of yi. We might as 
well consider the 2n values yi = Xi, yi = Xi + 1, 1 < i < n. By symmetry, it is 
enough to take i = 1. On the left, we apply Lemma [3.31 On the right, we observe 
that if 2/1 = Xi, all terms with 1 E S vanish, and if ?/i = Xi + 1, all terms with 
1 ^ S vanish. In both cases, the sum can be expressed as a sum over subsets of 
{2, . . . , n}. In this way, (15.11) is reduced to an equivalent identity with n replaced 
by n — 1. □ 

A particularly interesting case of (15.31) is 7 = \y \ — when it takes the form 
ry ( ^^ ^ [^ + \y\-\A+n-k] 

z^{x- 2/; A) = X. , 

with Ck independent of A. Moreover, ii r] = 1/N with N = 2, . . . ,n, we can even 
write (using that [x + N] = —[x]) 

[A + \y\ — \x\ + n — k\ 



Z^{x-y-X) = ^D, p^^^ 



The existence of such elliptic partial fraction expansions (in the sense of |Roj ) 
gives very precise information on Z„ as a function of A. 

Corollary 5.4. As a function of X, 

n 

^n(x;2/;A) J][A + j-l] 
i=i 

is a theta function of order n and norm \x\ — \y\ — (2)- Moreover, if rj = 1/N, 
with N = 2, . . . ,n, 

N-l 

Zn{x; A) JJ [A + n - j] 

j=i 

is a theta function of order N — 1 and norm \x\ — \y\ + n — (^) . In particular, in 
any case Zn has only single poles in X. 
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These facts are quite remarkable since, in general, the individual terms in (13 .Sp 
have poles of high multiplicity. The first part of Corollary 15.41 is also clear from 
Theorem 14.11 However, in the case t] = 1/N, the expression given there still has 
apparent multiple poles. 



6. Sums of determinants 

We proceed to show that in the important case = 1, the partition function 
can be expressed as a sum of — 1 determinants. 

When dealing with this case, we find it convenient to change to multiplicative 
notation. To this end, we observe that the expression 

is invariant under translations by l/rj in each of the variables Xj, t/i and A. This 
follows from Lemma [3l2] and Corollary 15.41 Thus, there exists a function Zn{x; y; A) 
such that 

...,q^";qy\...,qy-,q^)= q<\-\+\y\y^ Z^{x- y- A). 
We may write Theorem 14.11 as 



Zn[x] y; A) - 



^ y- TT diQy^ij)/yaii),y^ij)/xiq) -A- eiXq^'-^y^^j^/xj) ^^^^ 



Replacing also q'^ by 7 and recalling the notation (12.61) . Theorem 15. II can similarly 
be written 



Zn{x; y; A) 



(-i)(^)g(Ag") UlJ=lypi^^/yJ,1^^/yj) 

X V (-l)l^l^iMlS det (^^l^Hlll] (6 2) 



5C{l,...,n} e{Xq-\^\) l<^,J<n \ 9{xf/y,) 



where, in contrast to (15. 2p . 



X. 



s I Xiq, i E. S, 



Xi, i ^ S. 

We will need the following result, where we use the standard notation 

00 

{p;p)oo = Ylii-p'). 
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Lemma 6.1. When N is a positive integer, 



6{ax;p) _ {p^;p^)lo^i(^'':P) sr^ f^9{ax'^p^\p^^ 



9{x- p) {p; p)l d{x^- p^) ^ d{ap^- p^) 



Proof. Start from the Laurent expansion 



{p;p)'L0iax;p) 



OO f. 

Yl r^^' \p\ < 1^1 < 1' (6-3) 

^-^ 1 — an'' 



9{a,x;p) \ — ap 

which is a special case of Ramanujan's summation |GR^ Eq. (11.29)]. Replace 
the summation index /c by /c + A^j, where < < — 1 and j G Z. Observe that 
the sum in j can be evaluated using another instance of (16. 3p . The restriction 
\p\ < \x\ < 1 is removed by analytic continuation. □ 

We now assume that = 1, with = 2, 3, ... . Applying Lemma 16.11 with 
a = 7, X = Ag"~''^', the sum in (16.21) takes the form 



d{xl-i/yj] 



Ef-l 



X det 



i<ij<n \^ 9{xf/yj) 



X det /^ ^(^^^/^j) ^-k d{qixi/yj) 
i<i,j<n\ e{xi/yj) 9{qXi/yj) 



We thus arrive at the following result. 
Corollary 6.2. When q'^ = 1, 

{~lp{p^;p^)l9{\f 



Zn{x; y; A) 



(p; p)L e{q)^'e{^)--^e{x^; p^)Ye{xx^q^/Y) 

]\lj=iyp{xilyj,qxilyj) 

X - - 

l<'i<j<n 

„.^(7>^-'=-^ 

X 

A;=0 



ni<i<i<n Xjyjeixi/ Xj, yil yj) 
VAV'^^il^Vl^ Hpt ^^(WM _ .-^'^WM^ 



where 6{x) = 6{x;p). 

In particular, choosing 7 = p^^X~^ for some integer /c, one term in the sum 
vanishes and the partition function is expressed as a sum of A^ — 1 determinants. 
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For general g, the same method expresses the partition function as an in- 
finite sum of determinants. That is, if we use fl6.3p to expand the quotient 
^(A7g""''^')/^^(Ag"~''^') in ( 16. 2p we obtain the following identity. 

Corollary 6.3. Assume that \p\ < \Xq''\ < 1, < k < n. Then, 

(-l)(2)^(Ag") 



Zn{x] y; A) 



X 



ni<i<i<n XjyjO^Xi/ Xj, Vil Vj] 



oo 



7. Kuperberg's specialization 

Kuperberg showed how the partition function for the 6-vertex model can be 
specialized to the generating function (or t-enumeration) 

states 

where N denotes the number of entries equal to —1 in the corresponding alter- 
nating sign matrix. Using the Izergin-Korepin determinant formula, he computed 
the t-enumeration for t = 1,2, 3. We intend to generalize this for t = 1 and t = 2. 
To simplify the exposition, we first consider the case of general t. 

Lemma 7.1. Let, for each state, N denote the number of entries equal to —1 in 
the corresponding alternating sign matrix, that is, the number of blocks of the form 
i^l,--'). Then, 

Z^{q'^/\...,q''l^-l,...,l-X) 

0^)^g{3a+3b-c-d)/'i^ 



9{\q^) 



= g-KT)t-(^j) ^ n 

states blocks 

where the product is over all blocks ( ^ ^ ) and where 

t = 

^ e(gi/2)2- 

In particular, for p = 0, 
Zr.{q""\...,q-'"-l,...,l-X) 



1 - Ag" 

states blocks 
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where 

Proof. By definition, 

Z^{q-y^ . . . , 1, . . . , 1; A) = g^^^ H R^-ltliXn 

states blocks 

with local weights 

fl::(A) = fl::(A) = ,'/'^, 



We multiply all local weights by q~^^^0{q) /9(q^^'^), and accordingly the prefactor 
by (gi/^0(gi/2) /9(q))^\ We then muhiply Rlz by -g^/^ g^^^ _g-i/2, since 

these two types of blocks are equinumerous [Klj . this does not change the partition 
function. Using that there are factors of type and n + N factors of type 
R-+, we obtain 

Z^{q-'/',...,q-'/';l,...,l;X) 

\ y'il / statPS \ J / blocks 



where 



fl::(A) = -R::(A) = i, 



0(A) ' ' 0(A) ' 

Finally, one checks that in each case 



Q(^y^q{'ia+-ib-c-d)/A^^ 

0(Ag") ■ 

□ 
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8. Dynamical enumeration 

Kuperberg's proof of the alternating sign matrix theorem is based on the six- 
vertex model with q a primitive cubic root of unity. It seems interesting to consider 
the analogous specialization of the 8VS0S model. In the special case p = 0, we 
have the following result. 

Theorem 8.1. When p = and q = u = e^"/^, 
Zn{uj, . . . ,cc;; 1, . . . , 1; A) 

' (A„(l + ^"A2) + (-irc„a;2-A), (8.1) 



{1- Xu^+^){1- Xuj^+^j 
where 

(3j-2)! 



n + i - 1 ! 

is the number of alternating sign matrices of size n, and 

(3j-l)(3j-3)! 



Remark 8.2. C„ is the number of cyclically symmetric plane partitions that fit 
into a cube of size n. This was conjectured by Macdonald [M] and proved by 
Andrews [Xj. Moreover, the product AnCn equals the number of 2n x 2n half-turn 
symmetric alternating sign matrices. This was conjectured by Robbins [Ri| and 
proved by Kuperberg [K2]. 

Remark 8.3. If p 7^ and g = a;, it follows from Corollary 15.41 and generalities 
on theta functions that 



Zn{u), . . . ,u;; 1, . . . , 1; A) 

1 



e(Acu"+i, Ao;"+2) 

with Xn and Yn independent of A. We have not yet been able to find simple expres- 
sions for these functions, which seem to be natural elliptic analogues of An and C„. 
As will be explained below, a solution to this problem would compute the partition 
function for the three-colour model with domain wall boundary conditions. 

Proof of Theorem \8.1[ Let p = and = 3 in Corollary 16.21 Moreover, let 7 = 0, 
and replace each Xi by uiXi. The resulting identity can be written 



M{1 - ujY\1 - Acj"+i)(l - Acj"+2) 



Y\.l<i<j<rX-^0 ^i)iyj Vi 
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where 

Do= det ( i =^uj-u;'rx det ( 



= det — = (1 - ^ )" det 



l<i,j<n yi/j — UJXi Uj — Uj'^Xi J l<i,j<n yVj—X^ 

D,= det ( ^ e^')=(i_^)"y det 

l<i,i<n \yj — UJXi Uj — Xi J l<i,j<n — x^ 

It remains to let Xj, yi 1. This was done by Kuperberg for Dq (which is the 
same as D2), while Di can be treated by the same method. Indeed, it follows from 
|Klt Theorem 16] and also from [Ll Lemma 13] that 

lim — ; — det ' * 



x„...,x„^i Yli<i<j<ni^j - Xi)iyj - Vi) i^^'^^^" V y\ - 4 



It is easy to verify that the right-hand side equals 3 ^ ) when {k,l) = (1,3), 

and 3"^^ ^ when {k, I) = (2, 3). After simplification, using also that (1 — oo)^ = 
—3uj, we arrive at the desired identity. □ 

Consider now the case p = and rj = —2/3 of Lemma im so that g"^/^ = q = uj 
and t = 1. We find that the left-hand side of (18.11) equals 



UJ\ 2 ) 



1 - Xuj-^-'^ 



_ - Xuj'' 

states blocks 



Note that since d mod 3, the three numbers —c — d, c and d are all noncon- 
gruent mod 3, so 

1 - A^— ^ ~ 



Thus, we have 



EHtt 



(1-Ac<;«)(l-Aa;'=)(l-Aa;'^) 

states blocks ^ ' 

^ (1 - Acu")(A„(l + cu"A^) + (-l)"C„cu^"A) 
(1 - A3)"^+i 

Since each matrix entry a in the bulk of the height matrix gives rise to three 
factors 1/(1 — Acu"), the left-hand side equals 



5^ n (i-Acj'^)3 

states entries 
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times the correction 

(1 - x^Y-+yi - Xuj- 



(l-Acj2)2(i_Acj'^+i)2 

arising from boundary entries. Thus, we arrive at the following result, which 
clearly reduces to the alternating sign matrix theorem when A = 0. 

Corollary 8.4. For each state and for i = 0,1,2, let ki denote the number of 
entries of the height matrix congruent to i modulo 3. Then, 

Yi n (1 _ xu^f"^ 

states i=0 

_ (1 - Xu;')\l - A^"+^)^(A„(1 + iv-X') + i-irC^u'^X) 

This result is best understood in terms of 3-colourings. It was observed by 
Lenard jO] that reducing each entry in the height matrix modulo 3 gives a bijection 
from states to colourings of the (n + 1) x (n + l) square lattice with three colours, 
such that no adjacent squares have the same colour, and such that (in our case) the 
boundary condition arising from (13. ip is satisfied. It is then natural to introduce 
the partition function 



1 



2 • 



states 

In the case of periodic boundary conditions, the thermodynamical limit ?2 ^ oo of 
the corresponding partition function was computed by Baxter |Blj . Corollary 18.41 
evaluates the function for fixed n in the case of domain wall boundary conditions, 
when constrained by the relation 

1 1 lA^ 27 , , 

- + - + - = . 8.3 

Xq Xi X2 J XqXiX2 

Indeed, this surface is parametrized by = t/ (1 — Acj*)^, the dependence of t being 
trivial since kQ + ki + k2 = {n + If. The case of general p similarly corresponds 
to the unconstrained generating function. 

In spite of the constraint (18.31) . Corollary 18.41 contains much information about 
three-colourings. As an example, we pick out the coefficient of A on both sides of 
flO) to obtain 

{3ko + 3k^uj + 3k2UJ^) = -2A„(tu2 + tu"+^) + (-l)"C„w2'^. 

states 

Let 

states 

Since, for real values of a, b, c, 

a + buj + ccu^ = a = b = c, 
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we can deduce that 

3Ko - = 3Ki + 2An = 3^2 + 2A„, n = mod 3, 

3Ko = 3Ki = 3K2 + AAn - n = 1 mod 3, 

3A'o + 2A„ = 3iri-(-l)"C„ = 3K2 + 2A„, n = 2 mod 3. 

Together with the relation Kq + Ki + K2 = {n + l)^v4„, this allows us to solve for 
Ki. In Corollary 18.51 we state the result in terms of the probabilities 

(n + l)2A„- 

In other words, pi denotes the probability that a random square from a random 
3-colouring has colour i. We have also used that 

Cn ^ 2 ■ 5- ■ -(372-1) ^ (2/3)„ 
An l-4---(3n-2) (l/3)„' 

Corollary 8.5. J/n = mod 3, the probabilities pi are given by 



3 9{n + iy ' ' 9(n+ l)2(l/3)„' 

v =v - ^ ^ If ir+^ • 

3 9(n + l)2^^ ' 9(n + l)2(l/3)„^ 

if n = 1 mod 3, 

1 , 4 , , (2/3)„ 



Po = Pi = - + 777 — TT^ + (-1 



3 9(n+l)2 ^ ^ 9(72 + 1)2(1/3), 



3 9(n+l)2 ^ ^ 9(n + l)2(l/3)„' 
and if n = 2 mod 3, 



Po = P2 = 77 - 777 — + (-1 



.n+l (2/3), 



3 9(72 + 1)2 V ; 9(72 + 1)2(1/3), 



3 9(72 + 1)2 V J 9(72+ l)2(l/3)„' 

We point out that the asymptotics of pi in the thermodynamical limit 72 — 00 
are easily investigated using Stirling's approximation. In the first approximation, 
writing 

(2/3). _ r(l/3) ,3 ,3 
(i/3)„ - r(2/3) " ^ 
(note that the term of order r2~2/3 vanishes), we obtain in each case 

Pi = ^ + (-l)"Ci72-^/=^ + ^272-' + 0(72-=^) (72 ^ Oo), 

with explicit constants Ci and C2. 
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9. Dynamical 2-enumeration 

When q = —1, the partition function factors exphcitly. Indeed, both Theo- 
rem 14.11 and Corollary 16.21 reduce in this case to known elliptic determinant eval- 
uations. This can also be deduced from recent results of Foda et al. |FWZ] (as 
remarked in |FWZ^ §2.6] a special case of the Felderhof-type model considered 
there corresponds to the free fermion point, i.e. q = —1, of Baxter's model). 

Proposition 9.1. When q = —1, 

1 xea-ir+^XY/x) 



Znix;y;X) 



Yl Xiy^e{-Xj/x^,-yj/yi). 



X 

l<i<j<n 



First proof. We start from the case q = —1 of (16. ip . Since 

n n 

it can be written 



n 

T\e{-yi/xu)e{{-ir-^\y,/x,) n e{y,/xk) 

k=l k=i+l 

By a generalization of (12. 5p due to Tarasov and Varchenko |TV] . given more 
explicitly as [RoSl Corollary 4.5] (see also [HI Lemma 1]), the determinant equals 

n-l 

(_i)(2)x^-"e((-i)"+Ur/x) W 9{{-iy~^X) 

X W Xiyie{-Xj/x^,yj/yi). 

l<i<j<n 

Noting that 



9i-l) = 2- 



I oo 



2 

oo 



completes the proof. □ 

Second proof. Let = 2 in Corollary 16. 2[ and choose 7 = p/A^ so that the sum 
reduces to the term with = 0. The determinant is then 

f Ginxilyj) _ 9{-jXi/yj) 
i<i,j<n\ e{xi/yj) 6{-Xi/yj) 
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Applying the known identity 

6{ax) 6{—ax) 2x{p'^;p'^)l^9{a,apx'^;p'^) 

9{x) 9{—x) {p]p)1^9{x'^]p^) 
(which, for instance, follows from (16.31) ). the determinant is evaluated by (12. 5p . □ 

Combining Lemma 17.11 and Proposition 19.11 gives the following "dynamical 2- 
enumeration" . When A = 0, it reduces to 

5^2^ = 2©, 

states 

a result found already in jMRRj . 

Corollary 9.2. Let, for each state, N denote the number of blocks of the form 
( a-i "a^ ) • -^^'^ ^ = 0,1, 2, 3, let rrii denote the number of blocks ( ^ ^ ) with a + b + 
c + d = 2i (mod 8). Then, 



states 



2UJ 



r(i- A2)"'/2, 

(l+,A)(l-A2)("'-i)/2, 
[(l_,A)(l-A2)("'-i)/2, 



n 
n 
n 
n 



(mod 4), 
(mod 4), 
(mod 4), 
(mod 4). 



1 
2 
3 



(9.1) 



Proof. Let q = i in the second part of Lemma 17. H and evaluate the partition 
function using Proposition 19. 1[ This gives 



j^{c+d-3a-3b)/2^ 

1 - i-i^x 



2UV 



states blocks 

Since odd and even matrix entries interlace, 



1 + i"A 
1 + (-1)"A 



no 



-in) 



n even. 



blocks 



(l_A)«V2(i + A)"'/2, 
(l_A)("'+i)/2(i + A)("'-i)/2, nodd. 



independently of the state. Moreover, since a + 6 is odd, 

1 _ ^(c+d-3a-3fe)/2^ = 1 + ^{«+^+c+f^)/2 A 

This yields the desired expression. 



□ 



As an indication of the combinatorial meaning of Corollary 19. 2[ we identify the 
coefficient of A on both sides of (19.11) . We obtain 



2^(mo + imi — m2 — im-i) = 2(2) 



states 



0, n = (mod 4), 

1, n = 1 (mod 4), 
-2, n = 2 (mod 4), 
—i, n = 3 (mod 4), 
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or equivalently 



states 



'm2 - mo) 



2 -212;, n = 2 (mod 4), 
0, else, 



' 0, n even. 



states 
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